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OO ■ Abstract. If dark matter decays into electrons and positrons, it can affect Galactic 

radio emissions and the local cosmic ray fluxes. We propose a new, more general 
analysis of constraints on dark matter. The constraints can be obtained for any 
, decaying dark matter model by convolving the specific dark matter decay spectrum 

I with a response function. We derive this response function from full-sky radio surveys 

^ ' at 408 MHz, 1.42 GHz and 23 GHz, as well as from the positron flux recently reported 

Q^, by PAMELA. We discuss the influence of astrophysical uncertainties on the response 

^ ' function, such as from propagation and from the profiles of the dark matter and the 

. Galactic magnetic field. As an application, we find that some widely used dark matter 

C/3 decay scenarios can be ruled out under modest assumptions. 
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1. Introduction 

The existence of non-baryonic dark matter is supported by many astronomical and 
cosmological observations such as rotation curves of galaxies, gravitational lensing, 
cluster dynamics, large scale structure surveys and the cosmic microwave background 
(CMB) anisotropies [1]. However, the nature of dark matter remains elusive after 
decades of research - we only know that dark matter does not participate in the 
electromagnetic or strong interactions and that it has behaved as a non-relativistic 
fluid during the formation of the large scale structure of the universe. 

Many extensions of the Standard Model of particle physics generally predict new 
dynamics between the electroweak and the Planck scales together with a number of new 
particles, sometimes with the required properties to be dark matter. If dark matter 
is composed of these particles, it can decay on cosmological time-scales or, in the case 
that a conservation law forbids the decay, it can annihilate in pairs. Dark matter decays 
or annihilations in our galaxy will produce electromagnetic radiation, anti-matter and 
cosmic rays which we can be observed at Earth and serve as the first hints of dark 
matter non-gravitational interactions, provided their fluxes can be disentangled from 
astrophysical backgrounds. 

Such indirect detection of dark matter would not only provide us with a tangible 
glimpse of the existence of particle physics beyond the Standard Model, but would 
also allow a deeper understanding of structure formation. It may serve as a window 
into the early universe, probing times much closer to the Big Bang than primordial 
nucleosynthesis. In addition, indirect signatures of dark matter can reveal its large scale 
distribution and may thus serve as a diagnostic for astrophysics. Finally, decaying or 
annihilating dark matter can even have interesting consequences for stellar evolution. 

In this paper we develop a general formalism to derive information and constraints 
on decaying dark matter models. We focus on two observables, namely the cosmic 
positron flux on Earth and the synchrotron radiation from electrons and positrons in the 
galactic magnetic fields for both of which exist detailed observational data. The positron 
flux has recently been carefully measured by the PAMELA collaboration [21 [3] . For the 
synchrotron radiation, our analysis is based on three full-sky maps at 408 MHz [1], 1.420 
GHz [5] and 22 GHz [6] the latter from WMAP. 

Usually, constraints are derived in the literature for specific dark matter scenarios 
with given decay or annihilation spectra and branching ratios into final state 
products. For a given dark matter scenario one proceeds to compute the diffusion 
of electrons/positrons through the galaxy and the resulting synchrotron maps. One 
can then compare with observations and derive constraints as in [3 El El [ini [H] . This 
method can be tedious if one is to use detailed propagation models, specifically when 
including re-acceleration because a numerical treatment is then mandatory, or if one 
has to repeat many simulations to fit data or scan a parameter space. 

However, since the propagation equation is linear in the electron density, each 
injected electron energy evolves independently. Therefore, with a finite number of 
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numerical simulations at different injected energies we can construct a numerical 
response function of signal to background. These response functions do depend on 
astrophysical parameters such as the cosmic ray propagation model and the dark matter 
halo profile, but not on the microscopic decay scenario. Constraints can then be simply 
obtained by requiring that the convolution of the response functions with a given dark 
matter decay spectrum be smaller than unity. In this paper, we have developed our own 
numerical propagation code to build such response functions and provide constraints on 
decaying dark matter. As a service to the model-builder we also provide analytic fits to 
our numerical response functions. 

It should be noted that relatively large uncertainties from the propagation of cosmic 
rays and the dark matter halo density profile could make dark matter constraints 
relatively uncertain [121 HHj- We therefore systematically investigate the model 
dependence of our response functions by taking into account different choices for 
propagation models and dark matter density profiles as well as for the Galactic magnetic 
field. 

Recently, the PAMELA [2], [3] and ATIC [H] results show anomalously large fluxes of 
high energy positrons and electrons and indicate the presence of previously unaccounted 
nearby sources. The results from the PAMELA anti-matter search show a hard change 
of slope in the positron to electron ratio in the range 10-100 GeV while at the same time 
the observed anti-proton flux nicely satisfies the expectations from being a secondary 
product of interactions of cosmic rays with the intergalactic medium. The ATIC 
collaboration reported an interesting bump in the energy range from 300 to 800 GeV 
and a sharp cutoff near ~ 500 GeV in the flux of electrons plus positrons. However, 
more recent results from the FERMI/LAT experiment indicate at most a considerably 
smaller excess [15] than the suggested by ATIC. In order to explain these features, 
much discussion of the possible sources has ensued in the literature. Some conventional 
astrophysics sources such as a nearby pulsar [16], [TTJ [18], [JH] ED] can explain the positron 
"excess" rather naturally. However, given the impact a discovery of dark matter would 
have, much discussion has been on dark matter annihilation or decay as the origin of 
the PAMELA and ATIC "anomalies" , see for instance [2ll[22l[23l[2ll[25l[26l[271[28l[29] 
and [30] for a model independent discussion. Appealing as they are, these scenarios 
also face a number of potential challenges. For instance, most dark matter models 
generate a softer spectrum in electrons and positrons than the PAMELA and ATIC data 
imply [31]. In the dark matter annihilation case, boost factors of 10^ — 10^ are required 
to normalize the observed signals [21], but at the same time are in serious conflict 
with other observations such as radio emissions [32], 7— rays and anti-protons [23l [T2] . 
We believe that our study can help the model builder in the quest for explaining the 
PAMELA positron fluxes, or in providing complementary bounds from the synchrotron 
radio maps. 

The rest of this paper is organized as follows. In Sect. [21 we provide general 
estimates of positron and radio signals induced by dark matter, and describe our 
method of response functions. In Sect. [3] and Sect. [Uwe present the response functions 
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resulting from the observed radio emissions, and from the flux of electrons and positrons, 
respectively. In Sect. [5] we apply our method to some concrete dark matter decay 
scenarios. A summary and discussion are given in Sect. El In Appendix A we review 
the astrophysical parameters including propagation models, halo profile and magnetic 
fields and their influence on the distribution of electrons and positrons. In Appendix B 
we provide analytical fits to our response functions that can easily be applied to any 
dark matter decay model the reader may wish to probe. We will use natural units in 
which c = 1 throughout. 



2. Propagation of electron/positrons and Response Functions 

In this section we introduce response functions which describe the effects of propagation 
on the electron and positron fluxes and on the resulting synchrotron fluxes observable 
at Earth in terms of the injected electron or positron energy. To start, we give some 
simple estimates of these fluxes before embarking on a more detailed calculation of the 
response functions. 



2.1. Estimate of the electron-positron flux 

The energy loss time of electrons and positrons in a radiation field of energy density 
is 



-10"(^)s, (1) 

where the latter expression is often assumed throughout the Galaxy [33]. The 
confinement time due to diffusive propagation in the Galaxy is similar to the confinement 
time of hadronic cosmic rays which at GeV energies can be estimated from secondary 
beryllium isotopes in the Galactic cosmic ray flux [Mj . 

tconf(^) ^3 X 10^-10^^ s. (2) 

This is consistent with the diffusion time tdm{E) ~ /K{E) in a galactic disk of height 
2h ~ 4kpc with the diffusion constant [33] 

27 ^ \ „_2 „-l 



K{E) ^ 3 X 10^^ J cm^ s'^ , (3) 

which yields 

tME) - 3 X 10^^ (-^) ' s . (4) 



2kpcy VioGevy 

The effective lifetime of electrons and positrons is thus Te{E) ~ min [tioss(-E'), tdifr(-E')] 
At E ^ 10 GeV this is, therefore, re(10 GeV) ~ 10^^ s. 
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The differential ffux of electrons and positrons je{E) per energy interval from dark 
matter of mass mx and lifetime tx which produces on average Ye{E) electrons and 
positrons per decay and has a local density px can thus be estimated as 

E^UE)^E^^^r^iE) (5) 
Att mx Tx 

V0.3GeVcm-3; \ mx J \10^^ s J \ tx J 

Here, Ye{E){E/mx) < 1 depends on the particle physics of the decays and could be of 
order unity. 

The observed flux of electrons and positrons at ^ 10 GeV is [331 ES, fT5] 

E'^jf^E) ~ 2 X 10"^ GeV cm-^ sr"^ . (6) 

The gravitino dark matter model discussed in Ref. [36] was constructed such that 
the gravitino decays could explain the EGRET excess, leading to mx — 150 GeV, 
Tx — 10^^ s. Comparing Eq. ([5]) with Eq. ([6]) shows that the electron-positron flux 
produced by the decays can be comparable to or even exceed the locally observed 
electron-positron flux. A more detailed numerical simulation is, therefore, called for. 

Indirect effects such as radio emission in the Galactic magnetic field, can give 
complementary constraints since they are sensitive not only to the local electron-positron 
flux but also to the electron-positron flux induced by dark matter decay in remote parts 
of the Galaxy which is not directly measurable. We, therefore, now turn to radio 
signatures. 

2.2. Estimate of the radio flux 

The power per frequency interval emitted by an electron or positron of energy £" in a 
magnetic field 5, averaged over magnetic field directions, is given by 



P{u,E) = ^-^^^x' |i^4/3(x)Kv3(^) - liKj/six) - (7) 

where e and me are the electron charge and mass, respectively, and we have abbreviated 
X = z//(2i/c) where the critical frequency is 

.^(£,.±££(■^^ = 966^ f— ^VgHz. (8) 

' 47ime\mJ \6fiGj VlOOGeVy 

For u > 10 MHz self-absorption is negligible and the emitted radio intensity J(z^) in 
units of power per frequency interval along a given line of sight is then given by 

J(z/) = JdsJ dEUE)P{u,E), (9) 

where s is the distance along the line of sight. Inserting Eq. ([7]) with the approximation 
2x^{...} 6{2x -1.5) simplifies Eq. Q to 

o 7/2 1/2 r 
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where the critical energy 



47r 1/ \ / u \ 1/2 / ^ 



is the inversion of Eq. ([8]). Assuming the magnetic field approximately constant out to 
a distance d, for example for about 10 kpc towards the galactic anti-center, Eq. fllOl) 
can be quantified as 

2.6,10- {J^f (^) [JLy [E^UE)]U,, .(12) 
Inserting now the estimate Eq. ([5]), we obtain 

d \ f B 



/ 1/ \ 1/2 / 

.JM . 2.9 X 10-» (— ) ( 



lOkpcy \QfiG^ 
X , — ) (l^) (i^) ergcm-s-sr- 



(13) 



'Ye{E)E\ (r,{E)\ ^1026 s^ 

This is comparable to or higher than the measured high Galactic latitude radio fiux, 
which is of order 10~^ erg cm"^ s~^ sr~^ at GHz frequencies. We now turn to more 
detailed numerical calculations of the radio signatures and local fiux of positrons and 
electrons. 



2.3. Response Function and Constraints 

The propagation of positrons and/or electrons in the Galactic magnetic field is usually 
described by a diffusion model. Under this approximation, the diffusion-loss equation 
for the relevant particle density per unit of momentum interval n{r,p,t) can be written 
in the form 

On 

— -Vn = Q{v,p) (14) 
where the differential operator V is 



Vn = V- {D,,Vn - V,n) + ^ (p^D.^^^) ~^\pn- |(V ■ V,n) . (15) 



dp \^ ^^dpp'^J dp 

Here, D^x is the spatial diffusion coefficient, Vc is the convection velocity, re- 
acceleration is described as the diffusion in momentum space and is determined by the 
coefficient Dpp, and p = dp/dt is the momentum loss rate. Since we are interested in 
relativistic electrons and positrons, we will use energy E and momentum p indistinctly 
and write n(r, E) in the following. Henceforth, we use n+ and n_ for the positron and 
electron density, respectively, and = + n_. In the case of CP conserving decays 
one has n+ = n_ = ne/2. 

The source term for positrons/electrons due to decaying dark matter particles with 
mass and lifetime r is given by 

Q,(r,£„) = M£)^^ (16) 

mxTx dEo 
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where mx is the dark matter particle mass and tx its hfetime, pjs:(r) is the dark matter 
density profile in our Galaxy, and dN^/dE^ is the spectrum of positrons/electrons per 
dark matter particle decay. 

Consider the stationary solutions n^°{r,E) to the propagation equation for 
monochromatic injection of positrons or electrons at Eq, i.e. the Green's function 
satisfying 

- V nf> (r, E) = ^^^6(E - Eq) . (17) 
mxTx 

The solution of Eq. ( fT5l) for an arbitrary spectrum dN±/dEo can then be written as 

dN± , , 



n±(r,E) = y dE^nf{v,E)'- 



dEo ' ' 

The synchrotron flux arriving to the earth from a direction f2 characterized by 
galactic coordinates = (0, ff) has a contribution for each monochromatic injection E<^ 
given by Eq. ([9]) 

J^«(fi,z/) = ^jds J dEnf{Y,E)P{v,E) . (19) 

For an arbitrary injection spectrum dNe/dEo the synchrotron flux at frequency u 
is then obtained by 

J{Q,u) = j dEoJ^'{Q,u)^^. (20) 

It is then convenient to introduce the response functions for positrons Fp{E; Eq) and for 
synchrotron emission Fr{^l, u; Eq) as the ratio of the numerically computed n^\E) and 
J^°{(f), 9, u), respectively, and the observed fluxes as 

p.p. _ ng°(rearth,^) f \ f mx \ 

' °^ " nf^E) ll026sAl00GeVy 

= j^i^ (t5^) imkv) ^^^^ 

These functions depend neither on tx nor on mx and constraints on a given dark matter 
decay model can then be easily cast in the form 

/J....,.,..,^.(^)(^). 

The desired response functions can be computed numerically by using the methods 
of [391 ES]. In order to do so, we have developed our own numerical code. Details on 
our code and computations are described in Appendix A. 

Let us once more stress that our response functions do not depend on the specific 
decay spectrum, but still depend on the characteristics of the propagation model and the 
dark matter distribution. In this paper we use different halo models, see Appendix A, 
always normalized such that p(rearth) = 0.3 GeV cm~^. For other normalizations 
p(i'earth), our respousc functions have to be multiplied by p(rearth)/0.3 GeV cm~^. 
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Finally, apart from the injection energy Eq, the synchrotron response function 
Fr{fl, Eq) depends on the observed frequency v and the direction in the sky f2, whereas 
the positron response function Fr{E; Eq) depends on the observed positron energy E. 
In the latter case, we use the PAMELA data to get our constraints. It consists of seven 
different energies so we can construct seven different response functions, see Section 5. 
The synchrotron case is more complicated since in principle there are infinite directions 
to look at, and the optimal direction will depend on the injected spectrum dNf./dEQ and 
the observed frequency. We discuss this case in the next section. 

3. Response Functions for Radio Signals 

In this section we compute the radio emission induced by dark matter decay and 
establish the response function by comparison with radio observations. As can be seen 
from Eq. ([8]), the radio frequencies relevant for our study are between 0.1 and a few 100 
GHz. Although the cosmic microwave background (CMB) would dominate the radio 
sky at frequencies above ~ 1 GHz, this signal can be removed from the foreground 
based on the sensitive multi-frequency survey performed by the WMAP satellite. In 
Fig. [U, we show the full sky surveys at the frequencies 408 MHz [1], 1.42 GHz [5], and 
23 GHz [6]. We do not use the higher frequency channels (33 GHz, 41 GHz, 61 GHz, 
94 GHz) observed by WMAP, as they are considerably more noisy and less robust to 
foreground subtraction than the lower frequency bands. In addition, we smoothed all 
maps to angular resolution of 1°. 

3.1. Observed Radio Sky 

Based on various global multi-frequency model fits, recently a public code has been made 
available [ID] which allows to fit most of the radio survey observation in the range 10 
MHz-100 GHz to accuracy around 1% — 10% depending on frequency and sky region. To 
obtain the strongest possible constraints on radio emission due to dark matter decay, 
one would perform a pixel-by-pixel scan over the whole sky and over all frequencies 
between 0.1 MHz to 100 GHz. Most of the observed radio maps, however, have only 
partial sky coverage. Although these uncovered regions can be interpolated by global 
model fits, we want to rely on direct observations so we only use the 408 MHz, 1.42 
GHz and 23 GHz full sky survey maps for our study. 

Recently, an excess of microwave emission in the inner 20° around the center of our 
Galaxy has been revealed in WMAP bands between 22 and 93 GHz, after a subtraction of 
the free-free, dust and standard synchrotron emissions. This excess, dubbed as 'WMAP 
haze" [mils], is distributed with approximate radial symmetry and is uncorrelated with 
the known foreground^. The use of this haze map with subtracted "known" foregrounds 

X The origin of the haze is currently a hot topic of debate. In principle it could be explained by 
conventionally astrophysical sources such as pulsars [47] but an explanation in terms of dark matter 
annihilation has also been suggested [TOl HH |44] . 
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-3.37 ^^^^^^^^^^^^^^H^^ ^^^^^^^^^^H 



Figure 1. Maps of the radio sky at frequencies 408 MHz, 1.42 GHz, and 23 GHz, 
from the top left and moving clockwise [H [5l [6] . The color scaling is the logarithm to 
the base 10 of the flux in crg/s/cm^/sr. 



could further strengthen our constraints on decaying dark matter easily by an order of 
magnitude as has been shown already in the case of dark matter annihilation in |45j . 
We have, in fact, also performed this analysis and found the same order of magnitude 
improvement in the constraints than Ref . [l5] . In addition, the 22 GHz emission is likely 
dominated by spinning dust, as well as synchrotron and thermal dust emission. The 
synchrotron emission can be constrained from its polarized emission [l6]. However, it 
should be kept in mind that the astrophysical backgrounds themselves depend to some 
extent on the not very well known galactic magnetic field and cosmic ray propagation 
parameters. Therefore, in this paper we want to be conservative and stick to the 
real observations so we do not present response functions based on the haze maps or 
other background subtractions. The subtraction of foregrounds will be improved by 
forthcoming radio data from Planck and at low frequencies from LOFAR and, in a more 
distant future, from SKA. 

3.2. Radio emission from dark matter electrons 

In Fig. [2] we show the radio emission from dark matter decay for the five propagation 
models of Tab. lAll and for the three halo density profiles of Tab. IA2[ For the sake of 
illustration, we adopted a dark matter of 100 GeV, a lifetime of 10^^ s, and we use a 
decay spectrum dNg/dE = 6{E — mx) such that the total energy goes into one electron. 

From our numerical calculation illustrated in Fig. [2l it is clear that the largest 
uncertainty of synchrotron radiation comes from the propagation models. The average 
radio flux can differ by a factor of ten. For the MIN model, since the height of the 
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Figure 2. Model dependence of radio signatures at 1.42 GHz induced by decays of 
dark matter particles with mx = 100 GeV, tx = 10^^ s, for an injection spectrum 
dNe/dE = 5{E — mx)- Results for the five different diffusion models of Tab. lAll (from 
top to bottom: MIN, MED, MAX, DC and DR) and for the three dark matter halo 
profiles of Tab. R2l (from left to right: Kramers, Isothermal and NFW) are shown. The 
color scaling corresponds to the logarithm to the base 10 of the flux in erg/s/cm^/sr. 
Note that the color scale corresponds to the same flux range in all panels for convenient 
comparison. 



diffusion zone is smallest, most of the radio emissions occurs at low latitudes. In other 
propagation models the radio emission is more extended because of the larger diffusion 
coefficient and the larger scale height of the diffusion zone which leads to more dark 
matter decays contributing. 

Compared to the diffusion models MIN, MED and MAX, the DC and DR models 
always produce smaller signals over the whole diffusion zone. This is mostly due to 
the larger diffusion coefficient which allows the electrons to escape more easily from the 
diffusion zone corresponding to fewer confined electrons. Meanwhile, the power of re- 
acceleration, described by Dpp in Eq. (fTSil . is also weaker since it is inversely proportional 
to the diffusion constant {Dpp oc D~^), see Eq. f lA.2p . This implies that re-acceleration 
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plays an important role in propagation models. 

We also study the variation of the synchrotron emission due to different halo profiles. 
In general, decaying dark matter with the NFW profile produces the largest average 
diffuse radio signals due to the relatively steeper slope of the density distribution. For 
the other two profiles the emissions are comparable with each other. Since the decay 
rate is only proportional to the density, uncertainties from the halo profiles do not alter 
the resulting dark matter constraints significantly. We do not take into account any 
possible small-scale structure of dark matter halos since due to the linear scaling of 
injection rates with dark matter density it has much smaller influence on fluxes than in 
annihilation scenarios. 

3.3. Response Functions 

Assuming a propagation model and dark matter profile, the radio emission produced by 
dark matter decay can be obtained for any given decay spectrum. Then an excess map 
can be calculated in comparison with observed radio maps, defined as the map of the 
ratio of predicted to observed radio flux in a given direction. One can scan the whole 
excess map pixel by pixel until the largest excess is obtained. This pixel, therefore, 
corresponds to the optimal direction for observation. 

Figure [3] shows several examples of these excess maps at the frequencies 408 MHz, 
1.42 GHz, and 23 GHz, respectively. For the sake of illustration, we assumed dark 
matter particles with mx = 100 GeV with an NFW halo profile and decaying into one 
monochromatic electron or positron. The most important feature in Fig. [3] is that the 
best directions for dark matter constraints do not point towards the Galactic center 
region. Although the dark matter signal close to the center is always larger than 
elsewhere, the observed background flux overcompensates it. Moreover, the optimal 
direction is not the anti-center where backgrounds are the smallest. The optimal 
directions tend to be not far from the center, and most of them concentrate in the 
southern hemisphere as many complex components such as giant molecular clouds and 
the north polar spur inhabit the northern hemisphere. The location of the warmest 
color which indicates the largest excess not only depends on which propagation model 
and halo profile is assumed, but also depends on frequency. 

For constructing the response function, we have to perform different simulations 
with mono-energetic energy spectra at different injected energies Eq. For each of the 
resulting excess maps the optimal direction for observation is slightly different. We 
don't want to provide more than one response function per observed frequency so we 
have to fix one particular direction. In order to do so, we add up all the excess maps for 
different energies and search for the optimal direction. The selected direction is then 
optimized for a perfectly flat spectrum whereas it may be a bad choice for strongly 
peaked or hard spectra. Fortunately, our calculations show that different selections for 
the optimal directions do not change the response function dramatically, at most by a 
factor of two in the worst scenarios. 
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Figure 3. Excess maps, i.e. contours of predicted to observed radio flux, for 
decaying dark matter with mx ~ 100 GcV, tx — 10^^ s, and an injection spectrum 
dNe/dE = 5{E — mx)- Results for the five different diffusion models (from top to 
bottom) of Tab. and for three survey maps at 408 MHz, 1.42 GHz and 23 GHz 
(from left to right) are shown. Note the logarithmic color scaling for the excess, where 
warmer color indicates larger excesses. 



In Figure H] we present our response functions for the three synchrotron frequencies, 
showing the dependence on the propagation models and halo profiles. The optimal 
directions are shown in the caption. As indicated before, the uncertainty of the response 
function is dominated by the propagation model whereas the influence of the dark matter 
halo profile is small. This is because the best direction points far from the Galactic 
center, see Fig. El where the different halo models considered are similar. On general 
grounds, if the optimal direction is close to the Galactic center, the excess emission 
will be more sensitive to the halo profile since there the dark matter density is more 
uncertain in the absence of sufficiently high resolution numerical simulations. 
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Figure 4. The model dependence of the response function based on radio emission, 
F/, is shown. The response function based on the observed radio sky at 408 
MHz, 1.42 GHz and 23 GHz (from left to right), respectively, are given. The red, 
green, blue, magenta and black bands denote the MIN, MED, MAX, DC, and DR 
models of Tab. lAl) respectively. The width of the bands represents the variation 
within the Kra, Iso and NFW halo profiles of Tab. IA2I The optimal directions 
are (0,0) = (291°, -13.9°), (291°, -13.9°), (233°, 25°) for the three considered radio 
frequencies, respectively. Analytical fits to these curves are presented in appendix B. 

As shown in Fig. HI for the radio excess maps, the MAX propagation model always 
gives the strongest constraints. The DC and DR models, which exhibit similar behavior, 
one clearly sees an exponential cut off at low injection energies, whereas the response 
functions for the MIN, MED and MAX models are dropping more slowly with decreasing 
energy. This is not surprising since, in the latter case, re-acceleration shifts lower energy 
electrons to higher energies. The drop at low energies in these models is strongest 
at the highest frequencies at which re-acceleration of the corresponding higher energy 
electrons is less efficient. We note that in order to reproduce the observed B/C data, 
the re-acceleration zone in the MIN, MED and MAX models should be limited to a 
slab of height h^eac — 0.1 kpc, comparable to the height of the gaseous disk. If the re- 
acceleration region would extend to the full height L of the diffusive region, the response 
function would be flatter and its values would be higher by about a factor of 3 above a 
few tens of GeV. 

To illustrate these points we show in Fig. [5] electron spectra in the galactic disc 
at 1 kpc from the center and 0.2 kpc above the disk, for different propagation models 
and injection energies. When re-acceleration is included, we get a noticeable bump 
in the spectrum at a few GeV. Above these energies, the energy loss generated from 
inverse Compton scattering and synchrotron emission dominates over energy gain by re- 
acceleration, and below a few GeV, re-acceleration over compensates the energy losses. 
Thus a visible bump appears when the electrons accumulate in an energy region where 
re-acceleration and energy losses offset each other. It seems that the amplitude and the 
position of the bump is independent of the injection energy below a few GeV. The large 
amount of electrons and positrons accumulating in this bump region induce most of the 
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radio signals around GHz frequencies. That is why the shape of the response function 
drops more gradually in the lower energy region in the MIN, MED and MAX models 
compared to the exponential drop in the DC and DR models. 

In the DC model, due to the absence of re-acceleration, the energy spectra appear 
as a sawtooth shape as the number of propagated electrons above the injection energy 
have a sharp cutoff. Similarly, in the DR model, since the electrons and positrons can 
not gain enough energy from re-acceleration due to the larger diffusion parameter and 
the smaller Alfven speed, see Eq. ( ]A.2p . the propagated energy spectrum above the 
injection energy tends to zero rapidly. 

Another interesting property of the response function is that it tends to fall at high 
energies. This can be understood as follows: Higher energy electrons either loose energy 
more quickly, or, if their energy loss length is still larger than the half height of the 
diffusion zone, can propagate further and can thus escape from the diffusion zone more 
readily. The diffusion length can be estimated as Dxx{E)tioss{E) ~ a few kpc, and 
is comparable to the the thickness of the diffusion zone in the MIN and MED models. 
As a consequence, in Fig. [5] for the MIN and MED models, the propagated spectrum is 
indeed smaller by factors of a few for the highest injection energies. In the case of the 
MAX model, the half height of the diffusion zone is considerably larger than the typical 
diffusion scale so that the injection energy has a weak impact on electron spectrum 
below 100 GeV. 

We note that the choice of the diffusion models and the injection energy affect the 
shape of the electron spectrum significantly only below 100 GeV, where the influence 
of the diffusion mechanism on electron and positron propagation is still important. For 
energies above 100 GeV, energy loss dominates the electron spectrum. The spectrum 
in this energy range thus should only depend on the ISRF and magnetic field and not 
significantly on the diffusion parameters, as is confirmed by Fig. [51 The plots for the 
DC and DR models, for which re-acceleration is insignificant, confirm our qualitative 
analysis. For instance, for 10 TeV injection, a flat spectrum appears at energies above 
100 GeV, below which the spectrum steadily drops due to electron diffusion. 

The response function reaches its maximum around the critical energy of Eq. ffTTl) 
corresponding to the energy at which electrons emit photons of the considered frequency 
u (namely about 5 GeV for 408 MHz, 9 GeV for 1.42 GHz and 20 GeV for 23 GHz, for 
slowly varying magnetic field of ~ 5/iG strength in the regions of interest). 

We also note that the response function for the DC and DR models tend to give 
stronger constraints than the MIN scenario for the 23 GHz map. Since the thickness 
of the diffusion zone in the MIN model is only 1 kpc, there is no strong emission from 
directions far from the Galactic center, similarly to Fig. [2l As a result, the optimal 
direction at 23 GHz points to high latitude. 
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Figure 5. The electron or positron spectra versus energy at r = liipc, z = 0.2k;pc in 
the NFW halo model. The solid, dotted, short dash, long dash and dotted-short dash 
line represent an injection energy of 1 GeV, 10 GeV, 100 GeV, 1 TeV and 10 TeV, 
respectively. Color keys are as in Fig. U) 



4. Response Function for Positron Flux 

Recently, PAMELA reported a relatively large positron fraction in the electron/positron 
flux above 10 GeV [3]. Possible explanations include as yet unknown nearby 
astrophysical sources or the decay or annihilation of dark matter. However, decaying 
dark matter models can be constrained by requiring the predicted positron flux to be 
smaller than the observed one. This can again be expressed in terms of a response 
function along the lines of Sec. 12. 3[ In order to convert the positron fraction given 
by PAMELA data [3] into the positron flux, we multiply it with the latest e~^e~ flux 
observed by the Fermi Telescope [15] and note that the parametrizations for the Galactic 
electron flux in Ref. [331 ES] is larger than the new Fermi data by a factor of about 1.5 
below about 30 GeV. Compared with the PAMELA data, the statistical errors of the 
Fermi data is sub-dominant because of the flner energy binning and the smaller statistical 
error of the flux. Therefore, for the statistical error of the positron flux we take into 
account only the statistical error of the PAMELA positron fraction data. The Fermi 
data are well flt by a simple power law expression Jg = 172.37 s~^m~^sr~^GeV~^. 
The resulting "observed" positron flux is shown in Fig. [61 The strongest constraints 
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Figure 6. The positron flux observed at Earth as obtained by multiplying the e+ + e~ 
flux observed by FERMI [T3] with the positron fraction measured by PAMELA 
see text. 



come from the positron flux at the lowest energy where the statistical error is negligible. 

The response functions based on the positron flux measured by PAMELA [3] are 
shown in Fig. [71 We only consider the high energy region above 10 GeV where the solar 
wind has no signiflcant influence. Compared to the synchrotron response functions 
the positron response functions are generally larger and, therefore, in general lead to 
stronger constraints. 

We shall note however that this would not be the case in CP-non-symmetric dark 
matter decay models in which the positron flux is suppressed with respect to the electron 
flux. In this situation one should note that synchrotron constraints based only on the 
electron density are stronger than those based on the locally observed electron flux itself, 
since the electron flux is about ten times larger than the positron flux. 

Fig. [7| shows a prominent feature in the response function based on the PAMELA 
data in comparison with the response function based on the radio emission: It depends 
mostly on the diffusion model but little on the dark matter halo proflle. This is because 
high energy positrons mostly come from nearby sources within ~ 1 kpc where different 
halo proflles yield very similar dark matter densities. Further, the constraints are weakly 
affected by the magnetic fleld since energy losses are dominated by the background 
radiation flelds. 

In addition Fig. [7] shows that the response function cuts off below the energy at 
which the positron flux is observed in the DC and DR models where powerful re- 
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Figure 7. The response function Fp{E) resulting from the observed positron flux 
given by PAMELA for various energies at which the positron flux was observed. 
The model dependence is also shown. The color key is as in Fig. ID Analytical fits to 
these curves are presented in appendix B. 



acceleration is absent, as one would expect since electrons essentially can only loose 
energy in this situation. On the other hand, the response function tends to peak where 
the injection energy approaches the observed energy. Above that energy the response 
function gradually falls off due to the faster diffusion effects of higher energy positrons, 
similarly to the behavior of the radio based response function discussed before. The 
MAX scenario predicts the largest locally observed positron flux since the stronger re- 
acceleration in the MAX models shifts the predicted peak of the energy spectrum to 
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larger energies, similarly to Fig. O We should note that the amplitude of the bump in 
Fig. [5] in the MAX model is somewhat lower than in the MIN and MED scenarios due 
to the larger diffusion parameter. However, since this bump is shifted to larger energies 
in the MAX model, the MAX scenario still gives the strongest constraint. 

5. Constraints on Dark Matter Models 

The response functions developed in the previous sections can provide interesting 
constraints on any dark matter decay model. In this section we discuss two simple 
examples that can be applicable to a number of realistic scenarios. 

Let us first consider an extremely simple case, direct decay into and 
electron/positron pair X —>■ e~^e~ with branching ratio be- This will give conservative 
constraints for models that also allow decays into other fermion pairs such as [48l |49] . 
Moreover, it is certainly the simplest scenario that could account for the PAMELA 
excess [501 EO]- The spectrum of injected electrons is then dN±/dEQ = be6{EQ — rnxf^) 
and the constraints from the radio and positron fluxes are therefore 

T ^, /lOOGeV\^, T , /lOOGeV\^, 

These bounds are shown on the left panel of Fig. [H] for all the different propagation 
models considered. We see that in general the positron data provide stronger constraints 
on this particular model than radio data, unless the dark matter mass is smaller than 
the PAMELA energy window i.e. nix < H GeV, and at the same time there is no 
strong re-acceleration, as is the case in the DC and DR propagation models. 

As a more realistic case, we have considered the scenario developed in Ref. |49j . 
This framework involves an extra U(l) gauge symmetry under which the standard model 
fields have no charge, a hidden U(l). Despite being secluded, interactions of this hidden 
sector with the standard model particles are realized through a tiny kinetic mixing [5l] 
X with the hypercharge U(l)y. In particular, the hidden gauge boson A' couples to all 
hypercharged particles with an additional suppression factor x which can have extremely 
small values, cf. [521 [531 ISH [55]. If this hidden gauge boson accounts for the dark 
matter and has a mass around ~ 200 GeV it can decay into lepton, quark and W^-boson 
pairs. Furthermore, if the A' lifetime is of order 10^® sec, the 7— rays and positrons 
resulting from the decay can explain the EGRET and PAMELA excesses. Convolving 
our response functions with the positron spectrum provided in [19] we find that, based 
on the continuum component alone, this model can be ruled out in the MIN, MED and 
MAX propagation scenarios due to the strong re-acceleration effects. The direct decay 
into e~^e~ , has a branching ratio ranging from 0.05 to 0.14 in the mass range 100 GeV 
> niA' > 300 GeV and is directly constrained by the left panel of Fig. [^ . 

§ This scenario is somehow similar to the gravitino with R-parity violation described in |561l57j. whereas 
the decay branching ratio of hidden gauge boson to the W boson is highly suppressed compared to the 
gravitino case. 
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Figure 8. Bounds on a decaying dark matter particle for the decay mode X e^ + e~ 
(left panel) and X /+Z~+invisible (right panel), see the text for details. The bound is 
on the lifetime of the dark matter particle divided by the branching ratio of the relevant 
mode. The color key for the propagation model is as in Fig. [4l Constraints from radio 
emission are dashed and from the PAMELA positron flux (with normalization given 
by the newest FERMI data) are solid. Each constraint line is based on either the 
synchrotron response function for the three frequencies (408 MHz, 1.420 GHz and 23 
GHz) or on positron response function for the 7 PAMELA energy bins. 



Let US now consider the three body dark matter decay X ^ + /~+invisible 
where are Standard Model leptons and invisible stands for a nearly massless fermion 
{mx ^ 'Tiinv)- This decay arises for instance in supersymmetric versions of the hidden 
U(l) extension of the Standard Model outlined above, which also include the hidden 
gaugino, the supersymmetric partner of the hidden gauge boson |58]. As a concrete 
example consider the scenario recently described in [59]: Dark matter is made of bino- 
like neutralinos Xi which are the lightest supersymmetric particles of the Standard 
Model sector. The neutralino mass is taken to be around m^o = 300 GeV which gives 
the correct relic dark matter abundance while the hidden gaugino X is supposed to 
be a sub-dominant component of dark matter. If the hidden gaugino is lighter than 
the neutralino, the decay channel Xi ~^ l^l~X is available, whereas the reverse process 
X — > can be realized in the opposite case. These decays are mediated by heavy 

sleptons. Since squarks are usually heavier than sleptons, decays into quark pairs are 
suppressed. 

In order to study bounds on this decay we focus onto the democratic case, i.e. equal 
branching ratio for decaying into lepton pairs of the three different families. We used 
PYTHIA [60] to simulate the final electron and positron decay spectra, and our response 
functions to give the bounds shown in the right panel of Fig. [HI These decay channels can 
be also found in other contexts like in [61] where a massive B-L gauge boson mediates 
the kinetic mixing, in grand unified theories [62l [63], or in general SUSY models with 
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slightly broken R-parity with or without considering supergravity [MIESIES]- 
6. Conclusions 

In our study we have computed the prediction for two important signatures from 
decaying dark matter, namely the synchrotron emission in our Galaxy and the positron 
flux. In view of the recent experimental observations such as the radio full sky surveys 
and the new PAMELA data, we have introduced useful response functions that can be 
applied to constrain any interesting decay models. Robust constraints can be obtained in 
terms of convolving the response function with the specific decay spectrum into electrons 
and positrons. Our results show that the resulting constraints depend mostly on the set 
of propagation parameters rather than the halo profiles and reveal re-acceleration as an 
important ingredient. 

We have finally applied our method to provide model independent constraints on 
two widely discussed decay modes and studied the implications for some more concrete 
realizations. In some cases, specifically when choosing the propagation models with 
largest re-acceleration parameters our methods rule out recently proposed decaying dark 
matter models. 

We shall emphasize here that the methods showed in this paper can be used for 
other observables like gamma radiation, high energy neutrinos, etc... and can be also 
used for the case of annihilating dark matter. 

Dark matter decays or annihilations can also affect the extragalactic radio 
background. Compared with other astrophysical sources, dark matter could lead to 
discriminable anisotropic features in our radio sky [67]. Using future radio observation 
such as by SKA [68], possible signals of dark matter annihilation or decay could be 
discovered. In addition, new constraints could result from the 7— ray flux observations 
by the Fermi Telescope [69]. Moreover, decay would also affect the expansion history of 
the Universe because of the change of the equation of state [70], [71], and potentially leave 
an imprint on the Universe. Considering the heating and ionization effects on baryonic 
gas during the dark age [721 [73l [71] , future 21 cm observation [751 [761 [77] could discover 
visible evidence for dark matter decay or annihilations. 

Acknowledgements 

The authors want to warmly thank Christoph Weniger for fruitful conversations and his 
help with PYTHIA. This work was supported by the Deutsche Forschungsgemeinschaft 
(SFB 676 "Particles, Strings and the Early Universe: The Structure of Matter 
and Space-Time and GRK 602 "Future Developments in Particle Physics") and the 
European Union under the ILIAS project (contract No. RII3-CT-2004-506222). 

Appendix A. Numerical Solution of the Propagation Equation 
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In order to solve the diffusion equation for monochromatic injection of electrons, required 
for our response functions, we have developed our own numerical code. The main 
features are the same as in the GALPROP code [39]. We discretize the parameter 
space {v,p, t) using cylindrical coordinates for the position in the galaxy r = (r, z). The 
diffusion zone is confined to be a flat cylinder with radius rmax and height 2L. The 
number of bins used for the simulations in this paper is 60 in r and z and 80 in p (this 
last one in logarithmic scale). 

Neumann boundary conditions are imposed at the origin (r = 0, p = 0) since there 
is no net flux across these interfaces. In this point our code differs from GALPROP, 
since there no boundary conditions in p are used. Our code should then provide more 
accurate solutions at very low energies. 

We impose Dirichlet boundary conditions at the external surface of the diffusion 
zone by setting the electron/positron density to zero there. However, electrons are also 
produced outside the diffusion zone and some of them propagate into this region again. 
Because of this, the number of electrons in the stationary solution (within the diffusion 
zone) would be in reality a bit higher than our numerical results. In order to quantify 
this effect, a more consistent boundary condition is investigated Appendix A. 3. 

The stationary solution is looked for by using the Crank- Nicholson implicit updating 
scheme. The time intervals start with 10^ year and are decreased to refine the solution 
up to a minimum of 10^ years. We have cross-checked our results using GALPROP in 
a number of relevant examples and found the agreement satisfactory . 

In the rest of this appendix we provide details on the parameters entering the 
diffusion-convection-loss Eq. (fT5|) that we use to compute the response functions. 



Appendix A. 1. Diffusion, Convection, and Re- acceleration 

The diffusion term reflects the spatial propagation of cosmic rays through the tangled 
Galactic magnetic fields. The diffusion coefficient ^^.^.(r,^) is assumed to be constant 
within the slab considered and is described by using a rigidity dependent function, 

Dxx = PDo (^-^^ (A.l) 

where P = v/cis the velocity and R is the rigidity of the particle defined hjR = pc/Ze 
in terms of momentum p and electric charge Ze. The normalization Dq and the spectral 
index 6 can be determined from Boron-to-Carbon ratio data 



For case of re-acceleration the momentum diffusion coefficient Dpp is related to the 
spatial diffusion coefficient D^x using the formula given in Ref. [79] . 

where va is the Alfven speed, and w is the ratio of magnetohydrodynamic wave energy 
density to the magnetic field energy density, which characterizes the level of turbulence. 
We take w = 1 (since it can be subsumed in va)- The re-acceleration term Eq. (1A.2I) is 
restricted to a slab of scale height h^eac which is in general associated with the gaseous 



Galactic Signatures of Decaying Dark Matter 



22 



iVlOQcl 




n 


D 
11 


T 


T/ 


(IV d dZ 


T/ 


h 

"reac 






fknp2 /Mvrl 
[K.pL /iviyij 






K in / 1; 


Kin /t; / KTir* 


K rn /t; 


fknrl 


MIN 


0.85/0.85 


0.0016 


20 


1 


13.5 





22.4 


0.1 


MED 


0.70/0.70 


0.0112 


20 


4 


12 
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0.1 


MAX 


0.46/0.46 
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20 


15 


5 
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0.1 


DC 


0/0.55 


0.0829 


30 


4 





6 





4 


DR 


0.34/0.34 


0.1823 


30 


4 








32 
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Table Al. Typical combinations of diffusion parameters that are consistent with the 
B/C analysis. The first three propagation models correspond respectively to minimal, 
medium and maximal primary antiproton fluxes, abbreviated by MIN, MED. and 
MAX, respectively. In the DC model, the secondary e^, p and p fluxes flt the data 
well, and the DR model can easily reproduce the energy dependence of the B/C data. 



disk and, therefore, smaller than the scale height of the diffusive region [80], see Tab. lAll 
below. 

The convection velocity Vc is assumed to be cylindrically symmetric and to point 
in the z-direction perpendicular to the Galactic plane. The divergence of this velocity 
gives rise to an energy loss term connected with the adiabatic expansion of cosmic rays. 
The energy loss term p is due to interactions of the cosmic rays with the interstellar 
medium (ISM), interstellar radiation field (ISRF) and synchrotron radiation in the 
Galactic magnetic field. The ionization. Coulomb interactions, bremsstrahlung, and 
inverse Compton losses are also taken into account [39] and play an important role in 
case of re-acceleration. 

The largest uncertainties in the predicted fluxes come from poorly known 
propagation parameters, in particular the possibility of re-acceleration of produced 
electrons and positrons. The corresponding uncertainty can reach one order of 
magnitude. We list in Tab. lAll five different combinations of propagation parameters for 
the models MIN, MED, MAX, DC and DR proposed in Ref. [39l EH ESI EQ] , which are 
compatible with the observed B/C ratio. In the present paper, we will not study certain 
more extreme propagation models such as the ones discussed in Refs. [83l EH ES] which 
consider relatively large convection terms and anisotropic diffusion with coefficients that 
are different for the radial and the cylindrical directions. We leave the study of such 
models to future work. 

Appendix A. 2. Dark matter halo profile 

In the galactic center the radiation density is ~ lOeVcm"'^, thus the diffusion length 
with diffusion coefficient Eq. ([3]) during an energy loss time Eq. ([1]) is 

rfdiff ^ 0.45 (i?/GeV) kpc , (A.3) 

which at a distance of ~ 8.5 kpc from the Galactic center corresponds to an apparent 
angle 

Odiff ^3.1° (i?/GeV)-°-^ (A.4) 
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Table A2. Some of widely used dark matter halo density profiles used in the present 
study. 

This is smaller, but not much smaller than the off-set from the Galactic center of the 
typical sky directions we use to establish the constraint Eq. (l22l) from radio emission. 
We, therefore, expect a moderate dependence of the radio emission based constraints on 
any possible spikes concentrated at the Galactic center. On the other hand, dark matter 
profiles normalized to the dark matter density at the solar system differ significantly at 
~ 1 kpc from the Galactic center. Some dependence on the dark matter profile can thus 
not be avoided. 

In this work, we will adopt three spherically symmetric density profiles for the dark 
matter distribution. Recent N-body simulations suggest that a radially symmetric dark 
matter halo profile can be parametrized by [1] 

Po 



Px{r) 



(A.5) 



{r/roY'll + (r/ro)"](^-^)/" ' 

In Tab. IA21 we give the values of the parameters (a, /3, 7) for some of most widely used 
profiles such as the Kra [86], NFW [87] and Iso [88] profiles. The constant po can be 
normalized to the dark matter density of 0.3 GeV/cm^ at the local solar system (r = 8.5 
kpc). 



Appendix A. 3. Boundary conditions 

In solving Eq. f[T^ . traditionally, ones imposed the Dirichlet boundary condition 
ne{r,z = ±L) = 0, n^^r = r^^.^, z) = 0, at which the particles can freely escape. 
However, electrons and positrons are also produced by decays outside the diffusion zone 
in the Galactic halo. In this environment, they can propagate along straight lines since 
inverse ultra-relativistic inverse Compton scattering on low energy photons is boosted 
in the extreme forward direction. Furthermore, since ~ leVcm"^ and becomes 
eventually dominated by the CMB far above the Galactic plane, according to Eq. ([1]) 
the energy loss length is > 100 kpc up to TeV energies. Thus, energy loss can be 
neglected on halo scales. At cylindrical distance r from the Galactic center, the total 
flux from the halo into the diffusion zone at its boundary at 2; = ±L is given by 



je(E,p)halo 
f7r/2 



dN, 



X 



4:71 nixTx dE 

/•2tt 

d9 COS 9 sin 9 / d< 
Jo Jo 



{E) 

00 

ds px 



(A.6) 



r"^ + L"^ + P + 2s (r sin 9 sin (p + L cos 1 
where the dark matter profile pxi"^) is assumed to be spherically symmetric. 
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Figure Al. The dependence of electron density at Eartli on boundary conditions 
in tlie MIN propagation model, for injection of one electron of 500 GeV energy. The 
red line refers to the Robin boundary condition Eq. (|A.7|) . and the green one to the 
Dirichlet condition. 

In Eq. (lA.6p . the integration is performed over the hemisphere above or below the 
diffusion zone where the flux from a given direction is a line of sight integral over s. 
Continuity of the flux at the diffusive halo boundary then requires 

|D,,(E,r,±L)9,ne(^,r,±L)| = ^ne(E,r,±L)-j,(E,r)haio, (A.7) 

where ne{E,r,z) is the local electron plus positron density per unit energy whose 
distribution is simulated in the propagation code. Eq. flA.7p then determines the 
boundary condition for the numerical simulation of the electron-positron distribution. 
Although we use this third type boundary condition aX z = ±L instead of the Dirichlet 
boundary condition, the Dirichlet boundary condition turns out to be a very good 
approximation since most decays occur within 1 kpc from the Galactic center. For 
example, in the MIN model, considering the halo contribution, the radio signals would 
increase by only 10% in Fig. lAlt and in other propagation models with a larger diffusion 
zone the enhancement is negligible, as one would expect. 

Appendix A. 4- Magnetic Field and Photon Energy Density 

Various techniques have been applied to the determination of the Galactic magnetic 
field. Detailed analysis of the rotation measures and dispersion of pulsar emission has 
been carried out [89l [90] . The work presented in Ref . [91] which is based on the large- 
scale data set on starlight polarization [92] with nearly 7000 stars show that the local 
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Figure A2. The IRAS energy density as function of r at z = and as function of z 
at p = 0. Tlie contributions from stellar radiation, magnetic field from Eq. (|A.9p and 
from Eq. (jA.Sp as well as from the CMB are shown from top to bottom on the left 
side. 



field is parallel to the Galactic plane and follows the local spiral arms. A smooth 
Galactic magnetic field is also consistent with the conclusions of Ref. [91] and can be 
parametrized as 

5(r, z) = Q e(-"/20kpc)g-|.i/5kpc ^^_g^ 

Random fluctuations are not included in the model. In order to quantify the influence 
of uncertainties of the Galactic magnetic field on dark matter constraints, a second 
magnetic field model is also considered, which is parametrized by 

S(r, Z)=5 e(-('-8-5kpc)/10kpc)g-|.|/2kpc 

The value of these parameters are adjusted to match the 408 MHz synchrotron 
distribution [93] . 

The magnetic field profile close to the Galactic Centre is quite uncertain and could 
be considerably higher than a few fiG [9l]. However, this would further reduce the 
diffusion angle Eq. flA.4p and is thus unlikely to have a significant influence on the 
radio constraints since the discussion in Sect. 13.31 will show that the maximum excess of 
predicted over observed signal does not occur within ~ 5° of the Galactic center. 

The ISRF distribution can be derived based on the IRAS and COBE infra-red data 
as well as by using information on the stellar luminosity function. In the present work 
we use the model distributed with the GALPROP code [93]. In this model, the ISRF 
energy density is about 10 eV/cm^ near the center and 5 eV/cm^ at the solar position. 

In Fig. lA2l we show the spatial distribution of photon energy density including CMB 
and magnetic field. For the electrons and positrons with energy above 1 GeV, inverse 
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86 


10 


14 
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0.15 
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0.28 


0.29 


0.29 




6.59 


10.3 


31.2 


10.5 


12.4 


30 


4.5 


8.5 


31 


X 


2.27 


2.43 


3.94 


2.4 


2.5 


4.5 


2.6 


2.3 


3.7 



Table Bl. Values of the parameters in the fits Eq. (|B.f [) to the synchrotron response 
function for the different propagation models, assuming an NFW dark matter profile. 
The three frequencies ;/i.2,3 are 408 MHz, 1.42GHz and 23GHz, respectively. 



Compton scattering dominates electron energy loss since the stellar energy density u^-^d is 
always much larger than the magnetic field energy density ub- Since the synchrotron and 
inverse Compton processes have similar energy loss rates, the cooling time of electrons 
and positrons therefore are almost independent of the magnetic field. In other words, 
at very high energies where the diffusion length becomes large compared to the energy 
loss length, the number density of electrons and positrons is approximately determined 
by the strength of the interstellar radiation field only. 

Appendix B. Response Function Fits 

Appendix B.l. Response Function for Radio Signals 

The response functions for different halo models are very similar and in first 
approximation proportional to each other. For that reason, we will consider mostly 
the NFW halo in the following. The response functions for other halos are the same 
within roughly 10% . 

For the NFW profile the synchrotron response functions Fr{Eo) as a function of 
electron injection energy Eq in GeV can be fitted with the following expressions 

Fr{Eo) = Nc-^^'/^'^^'Eq^ (DC,DR) (B.l) 

Fr{Eo) = A^iE^^e-^-^*/^")" + A'2^o~^'(l - e-^^"/^*)") (MIN, MED, MAX) , 

where the fitting parameters are given in Tab. IBll 
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Table B2. Values of the parameters in the fits Eq. (jB.2p to the positron response 
function for the different propagation models and observation energies, assuming an 
NFW dark matter profile. 



DC 


10.17 


13.11 


17.52 


24.02 


35.01 


53.52 


82.55 


N, 


7.5 


9.1 


9.9 


8.9 


9.3 
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Appendix B.2. Response Function for Positron Flux 



The response functions for the positron flux based on the PAMELA data at a given 
energy as a function of injection energy Eq in GeV can be fitted with 



F°C(£;o) = N,{Eo/E - l)5ie-(^°/^) + N^E^'^ 

F?\E,) = ^^"^TsZ^,^ + (iVsi^o-'^e-^o/^ + N,E,^^) e^^^/^^)^ (B.2) 

where the values for the fitting parameters can be found in Tab. IB2I for the seven 
PAMELA energy bins E. 

The DC model does not include re-acceleration and can be modeled by a power law 
at high energies {Eq ^ E). Our calculations show an enhancement near E Eq which 
is fitted by the parameters with subscript 1. 

The DR model includes re-acceleration but the effect is not very strong. In the 
region Eq < E the response function is an exponentially suppressed power law while 
near the threshold and at high energies the behavior is similar to the DC model. 

The MIN, MED and MAX model have strong re-acceleration and can be fitted at 
both high and low energies with power laws with small indices. The matching at energy 
~ Eo is included by the exponential factors with an parameter x. 
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